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MET gODS AND SYSTEMS FOR ACCTMULATING 
HCS GENERATED BY A SEQUENCE ESTIMATION ALG( 



BACKGROUND 

5 The present invention pertains to methods and systems involved in sequence 

estimation algorithms and, more particularly, to methods and systems for efficiently 
accumulating metrics associated with a sequence estimation algorithm, e.g., as part of 
the process of detecting radio signals. 

The cellular telephone industry has made phenomenal strides in commercial 

10 operations in the United States as well as the rest of the world. Growth in major 
metropolitan areas has far exceeded expectations and is rapidly outstripping system 
capacity. If this trend continues, the effects of this industry's growth will soon reach 
even the smallest markets. Innovative solutions are required to meet these increasing 
capacity needs as well as maintain high quality service and avoid rising prices. 

15 In mobile communication, the transmitted signal is often subjected to a time 

smearing effect created by the time dispersive nature of the channel, i.e., the air 
interface between a base station and a mobile station. The channel effects are estimated 
in the receiver part of a communication system, and used by the detector to aid in 
attempting to correctly deduce die information symbols that were transmitted thereto. 

20 A commonly used technique for deducing such received information symbols is known 
as Maximum Likelihood Sequence Estimation (MLSE) which, implemented using the 
Viterbi algorithm, is optimal for situations involving Additive White Gaussian Noise 
(AWGN). 

An MLSE detector operates by selecting a known bit sequence closest to a 
25 received bit sequence. Because 2 k bit sequences are involved (k being the number of 
data bits within a frame) in a typical MLSE detector, the system stores the 2 k 
sequences for comparison with a received bit sequence. For a large value of k, this can 
be unwieldy . Viterbi simplified the maximum likelihood evaluation by noting that each 
of the states represented has only a finite number of possible predecessor states, and 
30 that only the path that agrees most with the received sequence (the minimum distance 
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path) need be retained for each state. Trellis diagrams are commonly used to illustrate 
the concepts of paths and states associated with MLSE detecting, an example of which 
is provided as Figure L For the simple four state trellis diagram depicted in Figure 1, 
the characteristic used by Viterbi to simplify maximum likelihood evaluation may be 
5 understood by considering that each of the four states (00, 01, 10 and 11) has only two 
possible predecessor states, e.g., you can only reach the 00 state from a previous 00 
state or a 10 state. Viterbi recognized that the likelihood evaluation could be 
significantly enhanced by using this factor in assessing possible paths through a trellis 
diagram. 

10 Each possible path, or combination of state transitions, has an associated 

likelihood, or probability, of being correct. The probability of any given transition is 
based on a newly received value, in view of a succession of predecessor values. These 
transition probabilities are commonly referred to as metrics, and a succession of 
metrics are referred to as a path metric denoting the likelihood of a sequence of 

15 possible state transitions. Exemplary algorithms for computing metrics are discussed in 
detail below. However, it will be appreciated by those skilled in the art that the 
accumulation of metrics during MLSE detecting is a rather computationally intensive 
function. Since processing resources continue to be a valuable commodity in today's 
electronic devices, system designers have sought for ways to decrease the 

20 computational intensity of MLSE detecting without sacrificing quality. 

One method for addressing this problenj is described in an article authored by 
G. Ungerboeck, entitled "Adaptive Maximum Likelihood Receiver for Carrier 
Modulated Data Transmission Systems," IEEE Trans, Commun., vol. COM-22, no. 4, 
pp. 624-535, May 1974, which applies two steps to reduce complexity. The first step 

25 is to expand the magnitude square term found in the metric computations and to 

eliminate terms that are common to ail hypotheses. As a simple example, the term (a- 
b) 2 can be expanded into a 2 -2ab+ b 2 . If "a" does not depend on the hypothesized data, 
then the a 2 term can be dropped from the metric computation. 

The second step described by Ungerboeck is to re-arrange the order of the 

30 metric computations. With standard MLSE detecting, metrics are computed and 
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updated based on successively received data samples. Each iteration of the Viterbi 
algorithm corresponds to a new received data sample. Using the second step as 
described by Ungerboeck, each iteration of the Viterbi algorithm corresponds to a 
newly transmitted symbol. 
5 Another approach to reducing the complexity associated with the MLSE 

detecting is described in U,S. Patent No. 5,499,272, entitled "Diversity Receiver for 
Signals with Multipath Time Dispersion", to G. Bottomley, the disclosure of which is 
incorporated here by reference. Therein, joint MLSE detecting and diversity 
combining can be performed using a variant Ungerboeck form by expanding the metric 

10 expressions and collecting terms that correspond to the same hypothesized symbol. 

Although MLSE detectors have been implemented through computationally 
efficient Viterbi schemes, e.g., using the Ungerboeck or variant Ungerboeck techniques 
described above, the MLSE detector might still be very computationally complex if 
large symbol alphabets are employed or if the number of required taps in the channel 

15 estimate is large. Both of these possibilities become more likely with next generation 
systems, e.g., EDGE which may employ 8-PSK modulation and a larger number of 
channel taps. 

Accordingly, it continues to be desirable to make sequence estimation 
calculations more computationally efficient and to design detectors that will have a 
20 reduced MIPS requirement so as to free computational resources for other 
functionality. 

SUMMARY 

These and other drawbacks and limitations of conventional detecting techniques 
25 are overcome according to exemplary embodiments of the present invention, wherein, 
according to a first exemplary aspect of the present invention, the addition of old 
survivor (i.e. , accumulated) metrics and delta metrics is reordered in a maimer which 
makes the calculations more computationally efficient. According to a second 
exemplary aspect of the present invention, symbol symmetries of modulation schemes, 
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e.g., phase shift keying (PSK) modulations, are used to make the calculations even 
more efficient. 

According to a first exemplary embodiment, the calculation of the delta metrics 
is partitioned into six metric terms (or four metric terms when the symbols have 
5 constant modulus), all depending on hypothesized sent data. As recognized by 

Applicants, one set of metric terms will be common for old states and another set of 
metric terms has no impact on the selection of survivor states. The former observation 
implies that one part of the delta metrics can directly be added to the old survivor 
metrics, while the latter observation implies that one part can be added to the survivor 

10 metrics after the selection of survivor states. In comparison with the known partial 
Ungerboeck approach, this exemplary embodiment significantly reduces the 
computational load associated with the accumulation of metrics, e.g., by approximately 
a factor of two when considering an exemplary radiocommunication system or terminal 
employing two tap channels and 8-PSK modulation. 

15 A generalized method for detecting according to the present invention includes 

the steps of: adding a first set of metric terms to a plurality of survivor metrics 
associated with a previous sample index to generate a first set of partially accumulated 
metrics; adding a second set of metric terms to each partially accumulated metric in 
said first set of partially accumulated metrics to generate a second set of partially 

20 accumulated metrics and selecting one of a plurality of states as a survivor state based 
on a comparison of said partially accumulated metrics in said second set; and adding a 
third set of metric terms to each partially accumulated metric in said second set of 
partially accumulated metrics to generate a plurality of survivor metrics associated with 
a current sample index. 

25 According to another exemplary embodiment, symmetries associated with the 

modulation employed can be used to reduce the amount of calculations even further and 
also decrease the storage requirements associated with metric accumulation. These 
effects are particularly significant when the hardware used to implement the radio 
stations, e.g., base stations or mobile stations, supports circular buffers and so-called 

30 butterfly instructions. 



WO 02/39689 PCT/SEO 1/02530 

-5- 

In existing products, the present invention will release computational resources 
previously employed for MLSE detecting that can be used for new functionality. In 
future products, the present invention implies that fewer, or cheaper, hardware units 
can be used and thus reduce the manufacturing cost. The present invention also 
5 permits unnecessary computations to be omitted such that other features, e.g., two 
pass equalization, can be performed resulting in large performance improvements. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The foregoing objects, features and advantages of the present invention, as well 
10 as other features, will be more readily understood upon reading the following detailed 
description in conjunction with the drawings in which: 

Figure 1 depicts a portion of a trellis diagram used to explain aspects of 
sequence estimation algorithms; 

Figure 2 shows a conventional radiocommunication system including plural base 
15 stations and a mobile telephone switching office; 

Figure 3 is a block diagram representation of a radiocommunication system 
which transmits and receives information over a time dispersive channel; 

Figures 4(a) and 4(b) depict exemplary state transitions and metrics associated 
therewith; 

20 Figure 5 is a flowchart depicting an exemplary method for accumulating metrics 

according to an exemplary embodiment of the present invention; and 

Figure 6 depicts oversampling associated with a fractionally spaced receiving 
design according to another exemplary embodiment of the present invention. 

25 DETAILED DESCRIPTION 

In the following description, for purposes of explanation and not limitation, 
specific details are set forth, such as particular circuits, circuit components, techniques, 
etc. in order to provide a thorough understanding of the present invention. However, it 
will be apparent to one skilled in the art that the present invention may be practiced in 
30 other embodiments that depart from these specific details. In other instances, detailed 
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descriptions of well-known methods, devices, and circuits are omitted so as not to 
obscure the description of the present invention. 

The exemplary radio communication systems discussed herein are described as 
using the time division multiple access (TDMA) protocol, in which communication 
5 between the base station and the mobile terminals is performed over a number of time 
slots. However, those skilled in the art will appreciate that the concepts disclosed 
herein find use in other protocols, including, but not limited to, frequency division 
multiple access (FDMA), code division multiple access (CDMA), or some hybrid of 
any of the above protocols. Likewise, some of the exemplary embodiments provide 

10 illustrative examples relating to the EDGE/GSM system, however, the techniques 
described herein are equally applicable to radio base stations and terminals in any 
system. Moreover, application of the present invention can be made to any 
implementation of sequence estimation techniques involving the accumulation of 
metrics, including those which are outside of the radiocommunication area. 

15 Figure 2 illustrates an example of a conventional cellular radio communication 

system 100. The radio communication system 100 includes a plurality of radio base 
stations 170a-n connected to a plurality of corresponding antennas 130a-n. The radio 
base stations 170a-n in conjunction with the antennas 130a-n communicate with a 
plurality of mobile terminals (e.g. terminals 120a, 120b and 120m) within a plurality of 

20 cells 1 lOa-n. Communication from a base station to a mobile terminal is referred to as 
the downlink, whereas communication from a mobile terminal to the base station is 
referred to as the uplink. 

The base stations are connected to a mobile telephone switching office (MSC) 
150. Among other tasks, the MSC coordinates the activities of the base stations, such 

25 as during the handoff of a mobile terminal from one cell to another. The MSC, in 
turn, can be connected to a public switched telephone network 160, which services 
various communication devices 180a, 180b and 180c. Both the mobile terminals 120a, 
120b, and 120m, and the base stations 170a-n can incorporate detecting structures and 
techniques according to the present invention. 
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Figure 3 illustrates schematically a radio transmission system for signal 
transmission over a time dispersive channel. A transmitter includes a unit 310 which 
receives an information carrying signal and generates corresponding digital symbols 
S(k). These symbols are delivered to a processing unit 311 which includes an output 
5 transmitter filter and a digital/analog converter. The symbols S(k) are signal processed 
in the unit 311 and output to a radio transmitter 312, which transmits the analog signal 
in the form of a signal R(T) at a selected carrier frequency. The analog signal is 
transmitted over a radio channel 313 to a radio receiver 314. The channel 313 subjects 
the signal R(T) to multipath propagation (among other things), as indicated in Figure 3 
10 by the double signal paths. The signals travelling along one signal path are reflected 
against, for instance, a building 319 prior to reaching the receiver. The radio receiver 

314 demodulates the received signal to baseband and delivers a baseband signal y(T) to 
a correlating and sampling circuit 315. In turn, the correlating and sampling 

circuit 315 delivers an observed sampled signal referenced y(t). The signal y(t) is 
15 received by a prefilter circuit 316 which outputs a prefiltered, observed signal z(k) to 
an MLSE detector 317. The signal z(k) is processed in the MLSE detector 317 in 
accordance with a sequence estimation algorithm, e.g., the viterbi algorithm, and 
delivers estimated symbols S D (/:), which should correspond to the transmitted symbols 
S(k). The MLSE detector 317 will determine the estimated symbols based, at least in 
20 part, on the accumulated metrics discussed below. The correlating and sampling circuit 

315 is also connected to a channel estimation circuit 318 and delivers thereto the initial 
values of a channel estimate which can be adaptively updated in detector 317 for use in 
the metric calculators. The selection of the correct proposed symbol sequence depends 
on the accumulated metric associated with each path tested by the detector. In this 

25 example, an accumulated path metric is generated for each symbol received. Each 
symbol received alters the assessment of the possible paths through the trellis. Hence, 
for each symbol received, a metric associated with the detector's level of certainty of 
the value, e.g., 0 or 1, of that symbol is conditionally added to the accumulated path 
metric. As mentioned above, it continues to be desirable to reduce the computational 

30 complexity associated with the task of calculating and accumulating metrics in MLSE 
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detector 317. To do so, it is useful to understand how these metrics are calculated at a 
more detailed level. 

To generalize, in a single branch case (i.e. no diversity) and with a mobile radio 
channel of L taps, defined by: 

h(0), h(l), h(L-l), (1) 

(where those skilled in the art will appreciate that the channel can be time-varying, 
although the time indices are omitted here to simplify the notation), the prediction error 
at sample index n of the sampled signal r(n) is formed as: 

e{n) = r(n)-^h{±)s(n-i) (2) 



where s(n-i) denotes hypothesized transmitted data associated with sample index n-4. 
10 With a symbol alphabet of length a, there will bea L prediction errors constructed from 

combinations of the sequence s(ri), s(n-l), s{n-L+ 1). The corresponding a L 

squared prediction errors are the metrics discussed above. 

Tfoie Viterbi algorithm (which is the exemplary sequence estimation algorithm 

used to illustrate the present invention in this specification) takes into account that the 
15 prediction error e(n-f-l) is independent of s(n-L+ 1). Hence, decisions associated with 

this hypothesis can be performed at sample index n. The a 1 " 1 states of the Viterbi 

trellis are formed from the combinations of the hypothesized data sequence inside the 

bracket in: 

2Q { s(n), .... s(n-L+2) , j(n-L+l)} (3) 

Into each state, there will be a paths of which only one will survive. Furthermore, 
each state k associates with an accumulated metric J n (k), k = 1, 2, a". 

Thus, at each iteration n % the following delta metrics (i.e., the metric to be 
added to the accumulated metric) are used in the Viterbi algorithm: 
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(4) 



where * denotes complex conjugation. By introducing the variables: 

g{i,i)=\hii)\ z t f(i)=r(n)h(i.)*, 1=0, (5) 
g(i t j)=h(±)hij)*, i=0, 1, 1-2, Jt-1, (6) 

the delta metrics (4) can equivalently be expressed as: 

jHr(n) | 2 +£gU, i) k(n-i) | 2 -2 £ Re{f (i) s (/i-i) *} 
+ 2£ £ Re{g(i,j-)s{ii-i) *s[n-j> } 



(7) 



1=0 



using the fact that x* y + xy* = 2Re{x* y}. As identified by Ungerboeck, and as 
5 applied in U.S. Patent No. 5,499,272, it is clear that the term \r(ri)\ 2 Ln equation (7) is 
independent of the symbol hypothesis. Thus, this term in equation (7) can be omitted 
to reach the simplified, "partial Ungerboeck" representation of the delta metrics. 

As can also be noticed from equation (7), the "partial Ungerboeck" 
representation can be further simplified for PSK modulation schemes where symbols 

10 have constant modulus, | s ( n -i ) | 2 =1 . Then the second term on the right hand side 
of equation (7) will be independent of the hypothesis, and can thus also be omitted. To 
better understand how metrics can be accumulated using the conventional partial 
Ungerboeck metrics, consider the following example for a communication system 
employing 8PSK modulation. 

15 Let M(k, p) denote the delta metric (the two first terms on the right hand side of 

(7) are omitted) associated with a state transition from state p to state k. The 
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accumulation of metrics for a conventional partial Ungerboeck scheme involving 8PSK 
modulation and a channel of two taps can be performed by the following pseudocode: 
for k = 1 to 8 
forp = 1 to 8 



+2 Re { g i 0 , 1 ) e u " lJ n/4 e n/4 ) 



(8) 



5 end 

J(*)=min{{J (1) 1) , . , . , J (8>+«(Jfc,8}>) ( q\ 

end 

However, Applicants have recognized that further simplifications are possible. These 
simplifications will now be described according to an exemplary embodiment of the 
10 present invention* Beginning with the partial Ungerboeck form (i.e., equation (7) with 
| r(n) | 2 omitted), the delta metrics are partitioned into a plurality of metric terms 
according to the present invention as: 

M / = A +A -2 (B +B +C +C ) CIO) 

12 12 12 

where: 

a - Lu,i)kM)| 2 (11) 

1 jl=0 



A 2 = g{h-X, L-l)|s(n-i+l)| 2 (12) 



15 
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r 

B = ^Re{f(i)s(n-i) *} (13) 
i i-o 

C = - £ £ Re{g(i,j)s(n~±)*s{n~j)} (14) 

1 i-Oj-i+l 



C = - 2^Re{g(i f L-l)s(/i-i) *sr(n-L+l> } (15) 

2 i-0 



Partitioning the metric calculations into these six terms permits them to be selectively 
calculated and added at different times relative to identifying survivor paths and making 
5 symbol (or bit) decisions as part of the sequence estimation algorithm, e.g., in the 

MLSE detector. For example, note that the terms A u B x and C x are independent of the 
hypothesis s(n~L+ 1), which implies that they have no impact on the selection of 
survivor paths. More specifically, A u B } and C x as partitioned above are scalar 
quantities that will be added to all of the possible paths and, therefore, only impact the 

10 size of the accumulated metric for the selected path rather than the identity of the path 
to be selected. Thus, Applicants have recognized that these terms can be added to the 
accumulated metrics after the selection of survivors;. 

Moreover, the terms A 2 and B 2 are independent of the states, which implies that 
they can be added directly to the old survivor metrics (i.e., the accumulated metrics up 

15 to sample index n-1), before state transition decisions are conducted. Partitioning the 
metric calculations as described above and recognizing these qualities of various metric 
terms, permits those terms to be recalculated less frequently than in conventional metric 
accumulation techniques. 

To further illustrate this aspect of the present invention, an example will now be 

20 provided which illustrates how the Viterbi scheme will operate when metric 

computations are partitioned and selectively accumulated as set forth above. For the 
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purpose of this example, consider 8-PSK modulation and L = 2, i.e. two tap channels. 
Those skilled in the art will appreciate that the present invention is equally applicable to 
different types of modulations and different numbers of channel taps. In this example, 
the number of states is then given by 8 W =8 and the total number of delta metrics to 
5 compute are 8 L =64. Thus, each state is associated with 8 delta metrics. For 8-PSK 
modulation, hypothesized data (symbols) belong to: 

s{n)e{e uk ~ 1)u/ \ k=l, 2, ..., 8} n=l, 2, W (16) 

A characteristic of PSK-modulation schemes is that the symbols have constant modulus 
| s | 2 = 1. This implies that the delta metric terms A , and A 2 are independent of the 
hypothesized data and can thus be omitted in the delta metric equation (10). This 
10 provides the following reduced expression for the delta metrics according to this 



example: 



M »= B +B +C + C 



(17) 



Note from equation (10) that minimizing the sum of M' 1 s is now equivalent to 
maximizing the sum of M //! s, with respect to hypothesized symbols for this 8-PSK 
example. 



15 



Given that A, and A 2 are not needed in this example, the partitioned delta metric 
computations for this example become: 



B x (k) = Re{f (0)e 



(18) 



B 2 (k) = Re{if(l)e 



-J t*-l)n/« 



(19) 



c (it) = o 



(20) 



C t {k f p) - -Re{g-(0, l)e 



(21) 
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where the state index p associates with the hypothesis s(n~l). This partitioning can be 
conceptually understood by referring to Figures 4(a) and 4(b). In Figure 4(a), it can be 
observed that the accumulated metrics for the state transitions have two terms in 
common, J„.i(l) and B 2 (l). Hence, these terms can be added before the total metrics 
5 are calculated since they appear for all transitions from state one to all other possible 
states. The same principle can be applied to state transitions from the other states. In 
Figure 4(b), it can be observed that the metric terms B,(l) and C x (l) are common for 
the states entering state 1. Hence, when comparing accumulated metrics, these two 
terms will have no impact at all on the outcome of the decisions, Le. , 

10 

(J (1) + B (1) +B(1) + C (1) + C (1,1)) 

fl-1 1 2 1 2 

-{J ( 2 ) + B ( 1 ) + B ( 2 ) + C ( 1 ) + C ( 1 , 2 ) ) 

n-1 1 2 1 2 (22) 

= J (1) +B (1) +C (1,1) -(J (2) +B C2) +C (1,2)) . 

«-l 2 2 n-1 2 2 

However, these two terms are part of the updated survivor metric J n (l), so they have to 
be added after the state transition decisions. Using these concepts, the accumulated 
metrics according to an exemplary embodiment of the present invention can be 
calculated according to the following pseudo code: 
15 1 . Before decision 

for p = 1 to 8 

^(P) ■ J a „ x (p)+B 2 (p) 

20 end 
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2. Decisions 
for k = 1 to 8 



J ( ic) = max ( { J (1)+C (Jt, 1), J (8)+C(Jr,8>) ^ 

end 

5 3. After decisions 

for k = 1 to 8 

* JJirJ+ajw+c^*) (25) 

end 

Above, the term is the accumulated metrics at sample index n-l. Since 
C 1 (fc)=0 it can of course be omitted in the pseudo code provided above, but has been 

10 included to emphasize its impact in general on the metric calculations. It will be 

apparent to those skilled in the art that by breaking up the metric accumulations in this 
manner, a significant reduction in the number of computations performed during the 
accumulation of metrics is achieved. 

From the foregoing pseudo code example, it will be appreciated that 

15 partitioning the delta metrics into the terms described above permits accumulating those 
metrics in various phases, which in turn results in a reduction of the computational 
complexity associated with MLSE detecting. A generalized method for detecting 
according to the present invention is illustrated by the flowchart of Figure 5. Therein, 
at step 410, only a first set of metric terms is accumulated before selecting a path 

20 through a particular state. In the foregoing example, this first set would consist of 
delta metric term B 2 , although for other types of modulation the first set may include 
other delta metric terms. 

At step 420, while decisions are being made, the method operates on a second 
set of metric terms, e.g., term C 2 . Lastly, at step 430, the method accumulates a third 

25 set of terms, e.g., J5, and C x (or just B x ), after path decisions are made. 
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When calculating all the delta metric terms in equations (11) to (16), the number 
of calculations can be further reduced by taking symmetries of the hypothesis into 
account. This will here be demonstrated for the 8-PSK modulation and two tap 
channels (L=2) studied in the above example. 
5 The 8-PSK symbol alphabet of equation (17) can be expressed as: 

. sell, je-' n/ \ j , -e-' n/ Vl, -ja- Wi , -J f e^'*} (26) 

Now, introduce the conjugated 8-PSK symbol alphabet 

s*e{l, e~ in/ \ j, -je~' n/4 ,-l, -e" in/4 , j, je" Jn/4 } (27) 

and the following %/4 de-rotated variables: 
TUT = f[i)e~ jn/ * t i = 0,l, g(0, 1) = g<0, 1) e~ jn/A (28) 

with^O and g(0, 1) being defined in equations (5) and (6). 

Next, consider B^k) and B 2 (k) specified in equations (19) and (20), respectively. 
10 By using (28) and (29), these variables can be expressed in a vector representation as: 



B ,u (1 > 










Re{jf(i) } 








TUT 




Re{TTIT} 


B ,o< 5) 










imifii) } 


= Re 




-J/tf) 




Xmtfli)) 
-Re{f(i)} 


B xn< 6 > 






-fU) 




-ReCTHT} 


B <n< 7 > 










-Imif(i) ) 


B ;n (8 > 






jTJJT 







Evidently, all the values readily follow from the real and the imaginary parts of 
the variables /(i) and 7(0 • If, for example, a digital signal processor (DSP) supports 
butterfly instructions, the number of loops (cycles) required in the calculations of J^fp) 



WO 02/39689 



PCT/SE01/02530 



-16- 

in equation (24), or J n (k) in equation (26), can be reduced with a factor of two, as 
illustrated below: 

f or p = 1 to 4 

J (p) ~ J (p) +B (p) 

end 

5 Continuing with delta metric term <7 2 (£, Z 7 ) specified by equation (22). By using 

equations (27) and (29), these variables can, in a vector representation, be expressed 1 



C t [k, 1) 
C 2 {k, 2) 
C 2 ik, 3) 
C 2 {k, 4) 

c 2 (k, 5) 

C 2 (k, 6) 
C (*, 7) 
C U, 8) 



from which we first note that: 

c z U, p) = ~c 2 (fc, p + 4) p = i/-»r4 (32) 



C{*,p> = -C 2 (lr + 4 # p) Jr«l, ...,4 (33) 

These symmetries assure that butterfly structures can be used in calculation ofJ n (k) 
10 in (25). By evaluating (32) for k = 1, it follows that: 







gr(0, 1) 








jfg(0, 1) 
jcrio, i) 




= -Re' 




-g(0, 1> 
-g(0, l) 


e -i(*-J|n/4 






-Jg(0, 1) 
-jgto, l) 








[ g(0,l>j 





(31) 
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C (1, 1) 




-Re{g(0, 1) > 


C (1, 2) 




lm(g{0, 1) } 


C (1, 3) 




lm{g(0, 1) } 


C (1, 4) 




Re{g(0, 1) } 


C (1, 5) 




Re{g(0, 1) } 


C (1. 6) 




-lra{g(0, 1)} 


C (1, 7) 




-lm{g(0, 1) ) 


C (1, 8) 




-Re{g(0, 1) ) 



(34) 



It can also be shown that: 

C 1) = C (q, 8) q = 1, 2, 3 

C <q+l, t+1) - C (g, t) t = 1, 2, ...,7 ( 35 > 

2 2 



This implies that if equation (35) is implemented in hardware as a circular buffer, all 
other terms p) can be detennined using different offsets rather than calculating 
each and every C 2 term directly. This will save both memory and computations, 
5 especially if the terms C 2 (k, p) would otherwise be calculated for every sample which 
will typically be the case when the mobile radio channel is considered to vary in time 
during a burst. 

While the foregoing examples are provided in terms of symbol spaced sampling 
and single antenna units, those skilled in the art will appreciate that the present 
10 invention can also be implemented in devices which include antenna diversity and/or 
fractional spaced receiver design. 

In fractional spaced systems and terminals, the sampling rate is a multiple U 
(the upsampling factor) of the symbol rate. For example, in Figure 6 a received signal 
r'(rc') has been sampled at the sampling instants n 1 - 1,2, with the multiple 
15 U = 2. By associating r x (n) = r'(2n-l) and r 2 {n) - r'(2n) a column vector of two 
elements can be formed according to: 

r(n) = \r x in) r (n) \ T n = l, 2,... (36) 
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When considering D antennas and fractional spaced sampling in general, the number of 
vector elements will be K = D x U. 

In the multi-branch case, the expression for the delta metrics of equation (4) is 
replaced by: 

AT (r(n) -£(n) ) *W(z{n) - *<ji) ) ^ 

5 where W is a KxK hermitian matrix that weights the different branches, r(n) is the 
received signal vector of K elements and t(n) > wftb elements defined by: 



Above, the * denotes both conjugation and transpose. 

By expressing equation (39) in a matrix form according to: 









f (n) 

K 




h (0) 

K 



= HS, 



(39) 



where the K x L matrix H contains the channel taps h m (i) and where the vector S 
10 contains the hypothesized data, the delta metrics in equation (38) can be written as: 

M = r*Wr + S*H*WHS-2Re{S*H*Wr] 



(40) 



By introducing the matrix: 



g(0,0) g(0,l) ... g(0,£-l) 
g(0,l>* g(l f l) 



gfO,L-l> * 



(41) 



and the vector: 



f(0) 
f(L-l) 



(42) 
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the delta metrics of equation (38) can equivalently be expressed as 

M~r*wr+S <7U, i) |s(n-i) | 2 -2 ^ Re [f { I) s (n-i) *} 

Jt-O 1-0 

(43) 

+2 S S Re{g(i f jf)*(ir-i) *s(n~j) } 

By comparing this equation (44) with equation (7), it is evident that they are of same 
structure. Hence, the same procedure partitioning and accumulation techniques 
described above for single branch, symbol spaced systems and terminals can also be 
5 applied to multi-branch, fractionally spaced terminals. 

In already existing products, the present invention will release computational 
resources previously employed for MLSE detecting that can be used for new 
functionality. The inventive techniques can be implemented in a processor, e.g., a 
DSP, associated with, for example, mobile terminals or base stations. In future 

10 products, the present invention implies that fewer, or cheaper, hardware units can be 
used and thus reduce the manufacturing cost. It also implies that unnecessary 
computations can be released such that e.g. , two pass equafetion can be performed, 
resulting in large performance improvements. 

The above-described exemplary embodiments are intended to be illustrative in 

15 all respects, rather than restrictive, of the present invention. Thus the present invention 
is capable of many variations in detailed implementation that can be derived from the 
description contained herein by a person skilled in the art. All such variations and 
modifications are considered to be within the scope and spirit of the present invention 
as defined by the following claims. 
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WHATIS CLAIMED IS: 

L A method for accumulating metrics associated with a sequence 
estimation process comprising the steps of: 

adding a first set of metric terms to a plurality of survivor metrics 
associated with a previous index to generate a first set of partially accumulated metrics; 

adding a second set of metric terms to each partially accumulated metric 
in said first set of partially accumulated metrics to generate a second set of partially 
accumulated metrics and selecting one of a plurality of states as a survivor state based 
on a comparison of said partially accumulated metrics in said second set; and 

adding a third set of metric terms to each partially accumulated metric in 
said second set of partially accumulated metrics to generate a plurality of survivor 
metrics associated with a current index. 

2. The method of claim 1, wherein said first, second and third set of metric 
terms are taken from the following delta metric equation: 

M=\r{n)\ 2 ^g[i, i) \s [n-i) | 2 -2 S Re{jf(i> s (n-i> *} 

I -0 I -0 

1-0 J-i+l 

3. The method of claim 2, wherein said first set of metric terms includes: 

A z = g(L-l, \s{n~L+l) | 2 

B 2 - Re{f[L-l)s(n~L+l) *) 
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4. The method of claim 3, wherein said method is employed in a 
radiocommunication system that uses a modulation having a constant modulus and said 
first set of metric terms includes only B 2 . 

5. The method of claim 2, wherein said second set of metric terms 
includes: 

C - - ^Re{g(i,L-l)s(n-i) *s(n-I,+l) } 

6. The method of claim 2, wherein said third set of metric terms includes: 

B = ^Re{f (i>s(n-i> *} 

1 1-0 



and 



c » - 



S -L Re{g(i, j)s(n-i) *s(n-j) } 

i-0 



7. The method of claim 6, wherein said method is employed in a 
radiocommunication system that uses a modulation having a constant modulus and said 
third set of metric terms includes only B x and C P 

8. The method of claim 1, wherein said process of sequence estimation is 
applied to detect symbols transmitted over a radio channel and wherein said previous 
index and said current index are sample indices. 

9. The method of claim 1 , wherein said sequence estimation algorithm is 
the Viterbi algorithm. 
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10. A method for processing a radio signal comprising the steps of: 

(a) receiving said signal on at least one antenna; 

(b) sampling said received signal to generate a sampled signal; 

(c) forming channel estimates from said sampled signal; 

(d) selecting at least one metric term based on a decision state of a 

detector: 

(e) accumulating said at least one selected metric term by adding at least 
one calculated value of said at least one metric term to a previously accumulated metric 
value; and 

(f) estimating a transmitted symbol sequence using a sequence 
estimation algorithm and said accumulated, selected metric terms. 

The method of claim 10, wherein said step of selecting further 

selecting a first metric term set before a decision is made by said 

The method of claim 10, wherein said step of selecting further 
selecting a second metric term set when decisions are being made by said 

The method of claim 10, wherein said step of selecting further 
selecting a third metric term set after decisions are being made by said 



11. 

comprises: 
detector. 
12. 

comprises: 
detector. 
13. 

comprises: 
detector. 



14. The method of claim 10, wherein said step of selecting further comprises 
selecting, as said at least one metric term, at least one of the terms from the following 
equation: 
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i-o l~o 

£ Re{g(i,j)s(n-i)*s(n-j)} 

i-0 



15. The method of claim 10, wherein said at least one antenna includes two 
or more antennas. 

16. The method of claim 10, wherein said step of sampling further 
comprises the step of: 

sampling said received signal with using a period that is less than a 
symbol period. 

17. The method of claim 10, wherein steps (d) and (e) are repeated for each 
of a plurality of decision states of said detector, with the selection of different metric 
terms in step (d) prior to reaching step (f), 

18. The method of claim 10, further comprising the steps of: 
calculating a first value of said selected metric term directly from an 

equation and, thereafter, determining subsequent values of said selected term by adding 
an offset to said first value, 

19. The method of claim 1, further comprising the step of: 
calculating a first value of each different metric term in said first set 

metric term directly from an equation and, thereafter, determining subsequent values of 
said selected term by adding an offset to said first value. 

20. The method of claim 1, further comprising the step of: 
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calculating a first value of each different metric term in said second set 
metric term directly from an equation and, thereafter, determining subsequent values of 
said selected term by adding an offset to said first value. 

21. The method of claim 1, further comprising the step of: 
calculating a first value of each different metric term in said third set 

metric term directly from an equation and, thereafter, determining subsequent values of 
said selected term by adding an offset to said first value. 

22. A device for estimating a sequence of values comprising: 

a processor which accumulates metrics as part of said process of 
estimating said sequence of values by performing the functions of: 

(a) adding a first set of metric terms to a plurality of survivor 
metrics associated with a previous index to generate a first set of partially accumulated 
metrics; 

(b) adding a second set of metric terms to each partially 
accumulated metric in said first set of partially accumulated metrics to generate a 
second set of partially accumulated metrics and selecting one of a plurality of states as 
a survivor state based on a comparison of said partially accumulated metrics in said 
second set; and 

(c) adding a third set of metric terms to each partially 
accumulated metric in said second set of partially accumulated metrics to generate a 
plurality of survivor metrics associated with a current index. 

23. The device of claim 22, wherein said first, second and third set of metric 
terms are taken from the following delta metric equation: 
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M=\r{n)\*+£g(±, i) \s (n-i) | 2 -2 £ Re{f (i) s (n-i) *} 
S Re(sr(i r j>5(i7-i) *s(iJ~7) } 



24. The device of claim 23, wherein said first set of metric terms includes: 

A 2 = g(l-l, L-l) |s(n-L+l) | 2 
B = ne{f(L-l)s(n-L+l) *} 

2 

25. The device of claim 24, wherein said method is employed in a 
radiocommunication system that uses a modulation having a constant modulus and said 
first set of metric terms includes only B 2 , 

26. The device of claim 23, wherein said second set of metric terms 
includes: 

C = - ^Re{g<i,L~l)s(n-i) *s(n-L+l) } 

2 1-0 



27. The device of claim 23, wherein said third set of metric terms includes: 
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1 i-o 



and 



c 

x 



)} 



28. The device of claim 27, wherein said method is employed in a 
radiocommunication system that uses a modulation having a constant modulus and said 
third set of metric terms includes only Bj and C v 

29. The device of claim 22, wherein said process of sequence estimation is 
applied to detect symbols transmitted over a radio channel and wherein said previous 
index and said current index are sample indices. 

30. The device of claim 22, wherein said sequence estimation algorithm is 
the Viterbi algorithm. 

31. The device of claim 22, wherein said device is a mobile terminal. 

32. The device of claim 22, wherein said device is a base station. 
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Fig. 2 
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Fig 4(a) 



1 0-=^ — — — — m-O 1 




J n - 1 (1)+B 1 (2)+B 2 (1)+ C 1 (2)-hC 2 (2,1) 



Fig 4(b) 



f Jn-1<1)+Bi(1)+B 2 (1)+C 1 (1)+C2(1,1) ^ ^ 




J n ^m+B t (1)+B 2 (V+Cf(1)+C 2 (1, 2) 
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Fig. 6 

) 

r'(n'-2) r'(n'-f) r'(n') r'(n'H) SAMPLE HATE 

~% • • : 

r t (n-1) r 2 (n~1) r f (n) r 2 (n) SYMBOL RATE 
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